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1. Introduction 

Since the emergence of Quantum Mechanics (QM), much effort was put into finding rela¬ 
tivistic wave equations that would govern the dynamics of mechanical systems. Schrodinger and 
then Klein and Gordon formulated a second-order wave equation, but at that time it seemed that 
the nature of the latter violated some fundamental properties of mechanical systems: the Klein- 
Gordon wave admitted a positive and a negative energy solution. British physicist Paul A.M. Dirac 
believed that the issue relied on the second-order nature of the differential equation. He therefore 
sought a first-order differential equation that respected the relativity principle. His theory was for¬ 
mulated in 1928 and the Dirac equation was later shown to describe relativistic spin 1/2 particles: 
fermions. This was followed by the development of Quantum Electrodynamics (QED), the rela¬ 
tivistic quantum theory of light and matter interactions which was then generalised into Yang-Mills 
(YM) theory, the theory of non-abelian SU (n) gauge fields. This summarises the success of parti¬ 
cle physics in the last century, success that culminated with the edification of the Standard Model 
(SM) of particle physics. Yet something can be seen as puzzling. Indeed, one of Dirac’s reasons to 
seek a first-order differential equation was the misinterpretation of the negative energy solutions. 
However immediately after the discovery of his equation, it became clear that Nature admitted par¬ 
ticles and antiparticles (positive and negative energy solutions). Nonetheless, fermions remained 
the only dynamical system that only admitted a first-order description. Indeed, both General Rel¬ 
ativity and Yang-Mills theory admit a first- and second-order formulation [3, 4]. As for fermions, 
because their equation is derived from the relativity principle, it also satisfies the Klein-Gordon 
equation; its Lagrangian is, however, first order in derivatives. Hence, it is natural to ask whether a 
(fully) second-order formulation of fermions that contains all the information of the original one is 
possible. 

Let us make the last argument more precise. A Dirac spinor can be written as the sum of two unitary 
infinite dimensional representations of the Lorentz group SO( 3,1) (or its double cover SL( 2, <C)): 

'Pa g (1/2,0) © (0,1/2) (1.1) 

which we call left-(unprimed) and right-handed (primed) respectively. The Dirac equation is then 
derived from the Dirac Lagrangian, here in 3 + 1 dimensions with the metric T] uv = (—, +, +, +): 

Sf D = W D (-i$-m) Vo (1.2) 

with ) = and we have the algebra of (Dirac) gamma matrices: 


{ydy v } = -2T^\ ( 7 A AW’, 75 ='WW, (75) + = 75 d-3) 

This Lagrangian generalises in a straightforward way so as to include the interaction of fermions 
and photons. We see that the Dirac equation 

(— i$ — m) 'P/j = 0 (1.4) 

relates spinors of one chirality to the other through the off-diagonal entries of the Dirac matri- 
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ces 1 . This is heuristically why a second order Lagrangian of the type 

&d = V D {-U + n?)'S D (1.5) 

does not work since the Klein-Gordon operator is diagonal and hence we lose information con¬ 
tained in the Dirac equation. 

A lot of insight can be gained on the issue when one expresses all the quantities in terms of 
two-component spinors (see for example [5]). In this paper, we will construct the simplest second- 
order spinor field theories; those of a (Weyl-)Majorana and of a Dirac fermion. Then, we will apply 
the same construction to the Standard Model. Finally, we will shortly discuss Perturbation Theory 
and Unitarity in this framework. 

2. Second-order (Weyl-)Majorana and Dirac fermions 

We start by considering the much simpler setup of a single (Weyl-)Majorana spinor and then 
repeat the procedure for a Dirac spinor. 

2.1 First-order Lagrangian 

The Lagrangian for a single Majorana fermion reads: 

Jz?Maj = —iV^X^dX — (, m/2)XX — (m/2)X^ X f (2.1) 

Flere Xa is a two-component spinor, X'' A is its Her mi t!an conjugate and 0 A ' A is the soldering form 

= t7,iiv (2.2) 

where r ]juv = diag(—1,1,1,1). We have also introduced the notation r) AA ' := (T' A/ ' <9 U and written 
the Lagrangian in an index-free way. Our index-free convention is that the unprimed fermions are 
always contracted as X% = X a %a, while the primed fermions are contracted in the opposite way 

W = W)a>(xY- 

2.2 Equations of motion and second-order theory 

The (first-order) equations of motion for the unprimed spinor are: 


—iV 2<9A — tnX' = 0 (2.3) 

The second-order formulation can then be obtained by carrying out the Berezin path integral 
over the primed spinors (which are treated as independent variables), which effectively amounts to 
substituting the above equation of motion in the Lagrangian. Then we have: 

Jz?Mai =- dXdX — (m/2)XX (2.4) 

m 

1 In the case of Majorana fermions the spinor is linked to its hermitian conjugate through the Dirac equation. 
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which is equivalent, after an appropriate rescaling of the fields to: 


j 2 ? Maj = -dXdX - {in' 1 /2)XX (2.5) 

The second-order Weyl theory is then simply obtained from the massless limit (after the path 
integral of the primed spinors has been solved): 


^Weyl = —dXdX 


( 2 . 6 ) 


2.3 Dirac fermions 

In order to describe Dirac fermions, we consider two massive Majorana fermions of equal 
mass. The system is then invariant under and internal SO(2) symmetry. Since SO(2) U(l), 
we can introduce a complex linear combination of the spinors and make the second symmetry 
explicit. It can then be made local by introducing a U( 1) gauge field and by minimally coupling 
the fermions. Thus, we define 

Dt; = (d-ieA)^ D X = (d + ieA) X , (2.7) 

where A^ = O^^'A^ is the electromagnetic potential and e is the electron charge. The gauge 
transformation rules are: for the fermions B, —> e‘ a q. X —> e~ ia X , and for the electromagnetic po¬ 
tential A^x —> An — (I /e)r)iiCX. The Lagrangian becomes 

•^Dirac = -iV2^D^-iV2 X t D X -m X ^-m^ X \ (2.8) 

where as before D := 6^D^. In order to obtain the second-order Lagrangian, the same procedure 
is applied to both primed spinors, and we obtain 


^Dirac = -2 D X DB, - m 2 X q , 


et 'V2 

=- D X , 

m 


X 


m 


(2.9) 


3. Short review of the Standard-Model 

For simplicity, we only consider the quark sector of the SM, for the complete picture, we refer 
the reader to [1]. 

3.1 Standard Model quarks 

The SM quarks can be put together in the following table 
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(1/2,0)reps 

SU{ 3) 

SU( 2) 

Y 

h 

Q = t 3 +y 

Qi = ( U> \ 

triplet 

doublet 

1/6 

1/2 

2/3 

w 

triplet 


1/6 

-1/2 

-1/3 

Ui 

triplet 

singlet 

-2/3 

0 

-2/3 

di 

triplet 

singlet 

1/3 

0 

1/3 


All quarks here are unprimed two-component spinors. Therefore, the Hermitian conjugate 
of Ui is denoted by u\, and Uj is another independent unprimed spinor. Notice that each SU(n) 
multiplet corresponds to a collection of unprimed spinors. Thus, e.g. Uj, which is a color triplet, 
has two types of indices suppressed: the usual spinor index, as well as the strong SU (3) index. 
With the index structure made explicit, this field is denoted by uf A , where A = 1,2 is the usual 
spinor index, and a = 1,2,3 is the index on which SU (3) acts. Only the flavour index / = 1,2,3 is 
left explicit in what follows. 

3.2 Higgs field 

This is the field that plays the central role in the the Standard Model. It is a complex scalar 
field of [7(1) hypercharge charge Y = 1/2. It is also a weak SU (2) doublet, i.e. it can be written as 
a column 

Higgs SU( 3) SU (2) Y r 3 Q = T 3 + Y 


Note that being an SU( 2) doublet, it is really a collection of 2 complex scalar fields </> ' and 
( p° (with Q = 0). We shall denote the weak SU (2) index by a,b,... = 1,2. Thus we can write the 
Higgs field as 0„, with 0i = 0 + and 02 = 0°- 


singlet doublet 


1/2 1/2 1 

1/2 - 1/2 0 


3.3 Quark sector of the Standard Model 

Using index-free notations, the Lagrangian for the quark sector of the Standard Model reads: 
quarks = — isflQ ' ‘DQj — isflii ' 1 Dili — is/ld "' 1 Ddj 

+ YWeQiUj - Yfa'Qidj - (y B + ) y fiJ Q)ef - (Y^dj ^ 

Here as before D AA ' = Qi lAA 'D^, where is the covariant derivative taht acts on each field 
according to its representation. The quantities Y‘ : are arbitrary complex 3x3 Yukawa matrices. 
The quantity £ = £ a b is the matrix 


= 



(3.2) 


which plays the role of a SL(2,C) metric. Then the object (j) T £Q = ($ T ) a eJ’ Qa is invariant 
under the action of the latter via Q —> gQ, <j) g(j) since g T £g = £. In particular, (j) T £Q is SU (2) 
invariant. It is then clear that all the interaction terms in (3.1) are SU (2) invariant. The hypercharge 
invariance is easily checked using the charges tables. 
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4. Second order formulation of the Standard Model 

We now formalise the procedure of integrating out the primed two-component spinors from the 
significantly more involved quarks Lagrangian. The equations of motion for the unprimed spinors 
are: 


Q] : isflDQ' = dj(yty/ 

u) : iV2Du‘= -(YjyjQ) (ep) 
dj : i\flDd l = -(yj)tfgj 0 


(4.1) 


Notice that some symmetry structure is making itself explicit in the equations of motion. Thus, 
let us combine the components of the Higgs field into the following 2x2 matrix: 


p^':=(afj) = 


( 0 °)* 0 + \ 

-<T 0 °) ' 


Under the weak SU (2) the matrix 4> transforms as: 

<f> f i-> £}<!>+, 


(4.2) 


(4.3) 


while the field p remains invariant. It is clear that p 2 = \(j)\ 2 is just the modulus squared of the 
Higgs field. 

In order to further simplify the equations of motion, a series of field redefinitions are needed, 
[1]. This leads to: 


Q] : iy/TDQi = -p & (&A). 
Q] : iy/lDQi = -p Q]& 


(4.4) 


Here, the doublet <2, transforms under the weak 5(7(2), and so does the Higgs matrix <f> : , 
while Qi does not transform. This suggests that we define new composite SU(2)-invariant quark 
variables <f><2, 


:= Q7- (4-5) 

This corresponds to a Higgs-field dependent SU (2) gauge transformation of the original quark dou¬ 
blet. As such, it can be pulled through the (gauge dependent) covariant derivative if one transforms 
the 5(7(2) gauge fields accordingly. The new vector field so defined will be an 5(7(2)-invariant 
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object. Notice that the SU (2) invariance does not disappear from the theory, it is simply “frozen” 
by the new set of variables that we chose as a convenient basis to work in. This change of variables 
for the gauge field is equivalent to working in the unitary gauge, but it should be stated that the 
second order formulation exists independently of the choice of variables. Keeping in mind this 
change in the derivative operator we can write the field equations as: 


Q] : iV2DQi = -p (e+A). 
Q] : isJlDQi = -p Q] 


(4.6) 


where from now on we drop the superscript inv from the Q, as it is clear that we will be working 
with this set of variables exclusively. We see that the equations become much simpler than in terms 
of the original fields. 

We now substitute the primed spinors obtained from the above field equations into the La- 
grangian (3.1) and obtain the following second-order Lagrangian: 

^quarks = - ^DQ'DQj - p (A Q) ' Qi, (4.7) 


where it should be remembered that the covariant derivative acting on Qi in (4.7) takes into 
account the field redefinition (4.5). 

Notice that this new second-order Lagrangian, even though it contains fewer terms that its 
first-order counterpart, is clearly non-polynomical in the (physical) Higgs scalar field p, because 
of the presence of 1/p in the kinetic term. In the case of Dirac theory (2.9) a constant rescaling 
of the spinors was all that we needed to obtain a canonical kinetic term. The same procedure can 
be applied to (4.7). However, p is now a dynamical field and thus, absorbing it into the fermion 
fields (again) changes the covariant derivative operators acting on both Qi,Qi. Denoting the new 
covariant derivative by St, we finally write: 


Quarks = -WQ^Qi ~ P 2 (AG)' Qi 


(4.8) 


where 1 / y ; p was absorbed into each spinor field. The new covariant derivative & contains non¬ 
polynomial Higgs-quarks interactions as well as the physical SU (2)-frozen gauge fields when act¬ 
ing on the unbarred doublet. It is clear that interaction vertices with the Higgs can be of arbitrarily 
high valency (due to non-polynomiality in p). 

The field equations (4.6) for the new fermionic fields of mass dimension one read 


iy/29Qi = ~P (2 + A)., iV2fft = ~pQl (4.9) 

and are now interpreted as the reality conditions, whose linearised versions are to be imposed on 
the external lines. 
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5. Perturbation theory and perturbative Unitarity 

We now discuss, in the framework of Dirac theory, how pertubation theory is modified but 
leads to the same results as in the first-order formalism, and also how Unitarity holds even though 
the Lagrangian is non-Hermitian. More details can be found in [2]. 


5.1 Feynman rules 

The propagator becomes very simple: 


{Q\T{Xa(p)£b{~p)}\0) =D(p) AB = 1 2 e AB 

p z + m 


(5.1) 


It is nothing but a scalar-type propagator, the espilon tensor in the numerator being the identity 
over the space of unprimed spinors. The complexity of the Dirac propagator is now shifted to the 
vertices. We have two interaction vertices with Feynman rules (incoming momenta): 


lie 


kc A £ba + Pb A £ca , 


—2ie 2 £ A B £Aii^ci) 


(5.2) 


Notice that the quartic vertex is simply the identity over both Minkowski spacetime and the space 
of unprimed spinors. 



Figure 1: Interaction vertices 


5.2 Spin averaged probabilities 

When we sum (or average) over photon polarisation states, one can make use of the Ward 
identities to obtain: 


— > Ptiv (5.3) 

pol. 


In our case, this will become: 


V , £aa ! £[] B i -a- — £ab£a'b' 

pol. 


(5-4) 
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As for the second-order fermions, we have: 


ejey ik)+e A el ( k ) 



(5.5) 


5.3 e fi —ten scattering 

Consider the simplest QED process: electron-muon scattering at tree level in the limit m e <C 
ntfi. Fig. 2. 



Figure 2: e jx —>• e p 


Let us first compute the amputated amplitude 'abcd for an incoming electron with momen¬ 
tum k\ scattered off an incoming muon with momentum p\. We have: 


^ABCD = 4 


(,ie) 2 (—i) 


q- 

4 ie 2 
q 2 




£de — hr/' £ae ) £ FF £e’F' \ P\b F £cf — Pic £bf 


(5.6) 


(&1 ■ Pi) AB £ CD — (&2 ■ P\) DB £-AC — {k\ • Pi) AC £ BD + {^2- Pi) cd £ AB 


where we defined: 


(k ■ p) AB := k A c 'p BC (5.7) 

and q 2 = (k\ —ki) 2 = (p i — pi) 2 = f- The complex conjugate amplitude is simply given by re¬ 
placing every unprimed spinor by a primed one and every primed spinor by an unprimed one, so 
that: 


Jt; 


4 ie 1 


A'B'C'D' 


(k i ' P\)a'B' £ c'd' — (^2 ■ P\)d'B' £ a'C ~ (k\ ■ Pi) A / C i £b'd' + ( ki • Pi) cd 1 £a'b' 

(5.8) 


and we defined: 


(k- p)A'B' '■ = kA'CPB C 


(5.9) 
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(k ■ p) AB e CD = m 2 k m 2 p 

(k-p)AB^CD < ( k-q) AC e BD = -\m 2 k {p-q) ( 5 . 13 ) 

(cj-1) cd £ ab =(k-p)(q-l ) 

where 

k- p = k A A ,p A ' A =k tl p B ( 5 . 14 ) 

Using this and neglecting terms proportional to the electron mass, we obtain: 

_ g^4 

\.£\ 2 = —r { h - pi)(k 2 -P2) + {k\-p2)(k 2 -p\) + m 2 (ki-k 2 ) ( 5 . 15 ) 

<r L J 

which is the well know squared amplitude for the unpolarised process. 

5.4 One-loop charge renormalisation 

Let us now look at a simple one-loop example. Although there arc half the number of fermions 
in our theory, fermions loops arc equivalent in both formalisms. Indeed, let us look at the amplitude 
for the charge renormalisation in the second-order formalism: 


10 



Second-Order Fermions and the Standard Model 


Johnny Espin 


iTT 


W 


A B 


= (- 1)4 e 2 j 


d 4 p 


P A 'b(p + k) B ' A + (p + k) A ' B p B ' A -\((p + k) 2 + p 2 ) £ AB '£ A b 


(2n) 4 

d 4 p £ a ' b 'eab 


[p 2 + m 2 } \{p + k) 2 + m 2 ] 


(5.16) 


+ (-l)4e 2 [ d p £ £AB 

’ J (2k) 4 [p 2 + m 2 } 

where we left explicit the contributions from both diagrams. 



Figure 3: One-loop photon two-point amplitude 


We can rearrange the latter to obtain: 

/n 1 -^ (k) A ' A B ' B 


= (- 1 ) 4 e 2 j 


d D p 


P A ' B(p + k) B ' A + (p + k) A ' B p B ' A+ m 2 £ A ' B ' Eab 


(5.17) 


(2k) 


D 


[p 2 4 - m 2 ] [(p + k) 2 + m 2 ] 


which is the usual tow-component version of the loop integral that has to be calculated in the 
first-order formalism. 


5.5 Remarks on perturbative Unitarity 

In quantum field theory, the S-matrix is generally written in the interaction picture: 


S = &j!* xSe *V> (5.18) 

where S denotes time-ordering. Taking matrix elements between physical states, Unitarity of 
the S-matrix reads: 


(/!*')= L (f\S\n)(n\sYi) = £ (/l 5l »H 5 1*'> 

phys n phys n 


(5.19) 
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where the sum runs only over physical intermediate states, which arc eigenstates of the free 
Hamiltonian, i.e. they correspond to on-shell particles. For perturbative computations, one splits 
the S-matrix into the identity plus a transition matrix T : 


S= 1 +/T 


(5.20) 


The Unitarity equation then reads: 


(f\(iT)\i) + (f\(iT)'\i) = - £ (f\(iTY\n)(n\(iT)\i) (5.21) 

phys n 

This equation is generally proven in Perturbation Theory, order-by-order in the coupling con¬ 
stant appearing in the interactions, and diagram-by-diagram. Usually, this is shown to be true as a 
consequence of the Hermicity of the Lagrangian, however, in [2], we show that this latter require¬ 
ment is not necessary and we prove Unitarity of second-order QED at all orders. The proof relies 
on the existence of a non-trivial real-structure that involves a derivative operator: 



This is nothing but the first-order Dirac equation imposed on the physical states of the theory. 
Imposed linearly (without any coupling to the gauge fields) on the free theory, this leads to a 
positive definite Hamiltonian Hq ~ a'a. While, imposed on the external states, they lead to a 
Unitary S-Matrix. 

6. Conclusion 

In this short note, we demonstrate that a second-order description of spinor fields is viable and 
that in many aspects, it is simpler than the usual first-order description. This new reformulation of 
the SM leads in particular to new insights on the available Unification patterns that could be allowed 
[1]. Furthermore, since the formalism involves only half the number of spinor fields and because 
the interaction vertices seen as matrices live in a much smaller space than their four-dimensional 
counterparts, the complexity of perturbative calculations is greatly decreased. It is also believed 
that for the purpose of lattice calculations, working with a scalar-type propagator could simplify 
the implementation of the theory. 

Finally, it should be noted that the construction that is followed here intends to mimic the SM 
calculations, but it might be possible to take the second-order formulation as fundamental and try 
to implement modifications of the theory in the latter. 
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